Relationship between spin squeezing and single-particle coherence in two-component 
Bose-Einstein condensates with Josephson coupling 
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We investigate spin squeezing of a two-mode boson system with a Josephson coupling. An exact 
relation between the squeezing and the single-particle coherence at the maximal-squeezing time is 
discovered, which provides a more direct way to measure the squeezing by readout the coherence 
in atomic interference experiments. We prove explicitly that the strongest squeezing is along the 
J z axis, indicating the appearance of atom number-squeezed state. Power laws of the strongest 
squeezing and the optimal coupling with particle number N are obtained based upon a wide range 
of numerical simulations. 
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I. INTRODUCTION 

Spin squeezing is a nonclassical effect of collective 
spin systems [H, El, E| 0], showing reduced spin fluc- 
tuation in one certain spin component normal to the 
mean spin. Kitagawa and Ueda proposed spin squeez- 
ing generated by the self-interaction Hamiltonian H± = 
2k J 2 , due to the so-called one-axis twisting (OAT) effect 
Q. The OAT-type spin squeezing could be realized in 
weakly interacting Bose-Einstein Condensate (BEC) [5|, 
or atomic ensemble in a dispersive regime @. The self- 
interaction Hi also leads to phase diffusion of the BEC 
7 1 .which indicates a decay of single-particle coherence 

1, 1,0, [H El O- 

Beyond the OAT model, an Josephson-likc coupling 
(JLC) term f2 J x was added to the Hamiltonian H\ with 
purpose to coherently control the phase diffusion [l3| and 
the spin squeezing [3, [3 EH ■ It was shown that the 
JLC model [see Eq.(Q}] results in strong reduction of spin 
fluctuation along the z (i.e., J z ) axis, provided that the 
additional field is tunned optimally [15| . We found the 
maximal-squeezing time to of the JLC model, and pro- 
posed a simple scheme to store the strongest squeezing 
along the z axis for a long time [16j |. So far, there remain 
certain questions unsolved: Is there any relation between 
the squeezing and the single-particle coherence? In addi- 
tion, to what degree can the strongest squeezing reach in 
the JLC model? The first question is important because 
it relates to measurement of the squeezing. 

In this paper, we present an exact relation between 
the squeezing and the coherence by solving the Heisen- 
berg equation. Our results show that local minima of 
the squeezing and the coherence occur simultaneously for 
the coupling f2 larger than its optimal value f^o- Un- 
like the OAT scheme, where number variance AJ Z is 
time-independent, we prove explicitly that the squeez- 
ing at time to is along the z axis in the JLC model 
171 ]. The strongest squeezing obeys the power law 
io = AJ z (t )/^/J/2 cx N- 1 / 3 , which can be measured 
by readout the single-particle coherence through the vis- 



ibility of the interference fringe @, H, El El, El ■ 

Our paper is organized as follows. In Sec. II, we in- 
troduce theoretical model and derive some formulas for 
the single-particle coherence and the squeezing param- 
eter. In Sec. Ill, quantum dynamics of the coherence 
and the squeezing are investigated for the OAT and the 
JLC models, respectively. In Sec. IV, we present exact 
relation between the coherence and the spin squeezing at 
the maximal-squeezing time to- In Sec. V, power rules 
of the optimal coupling and the strongest squeezing as 
a function of particle number TV are investigated based 
upon a wide range of numerical simulations. Moreover, 
we compare numerical result of to for the optimal cou- 
pling case with its analytic solution. Finally, a summary 
of our paper is presented. 



II. THEORETICAL MODEL AND SOME 
FORMULAS 

To begin with, we consider a two-component weakly- 
interacting BEC consisting of 2j atoms in two hyperfine 
states |1) and 12) coupled by a radio- frequency (or mi- 
crowave field lIlGl. A tightly-confined BEC can be 
described by the JLC Hamiltonian (h = 1) [2(| 



H 2 = SIJ X + 2kJ 7 



(1) 



where the angular momentum operators J+ = (J-y = 
a\a\, J z = {a\a,2 — a\a\)/2 obey the SU(2) Lie algebra. 
The total particle number N = a[ai+a\a,2 is a conserved 
quantity. In Eq.([T]), we have neglected the term propor- 
tional to J z by assuming equal intraspecies atom-atom 
interaction strengthes [j| . The Rabi frequency O can be 
controlled by the strength of the external field. The self- 
interaction term 2kJ z leads to spin squeezing, which is 
quantified by a parameter [l|: 



,•1/2 



(2) 



2 



where j = N/2 and (AJ n ) m ; n represents the minimal 
variance of a spin component J n = J ■ n normal to the 
mean spin (J). The coherent spin state (CSS), defined 
formally as [2l| 



Pm'(0) = Pm (t) = 0, i.e., c_ m (f) = -c m (t). Quantum 
dynamics of the odd iV system depends on the equations 
of pm . The above processes have certain advantages to: 
(i) reduce the total Hilbert space dimension from 2j + l to 
j + 1 (even AT) or j + 1/2 (odd N); (ii) Since c_ m = ±c m , 
we obtain (J a ) = (J 2 ) = and (J x ) = (J + ) 7^ 0, i.e., the 
mean spin is always along the x axis. Actually ( J+) is a 
real function; (iii) The correlation function is simplified 
as = j -1 1(^)1: and the spin component normal to 
the mean spin is J n = J • n = J^sin^ — J z cos9. The 
variance of J n is (AJ n ) 2 = (J^)-(J„) 2 = (C*- A cos 20- 
Bsin20)/2, where A = (J 2 — J 2 ), B = (J 2 J a + J y J x ), 
and C = (J 2 + J 2 ). From the relation <9 e (AJ n ) 2 | e = 0, 
we get tan(20 m ; n ) = B/A and the minimal variance 
= (C - VA 2 + B 2 ) /2. 



_ e —i&(J x sin<; 



- Jy COS E 



(3) 



has the minimal variance (AJ n ) m i„ = Wj/2 and £ = 1. 
Therefore, a state is called spin squeezed state if its vari- 
ance is smaller than that of the CSS, i.e. £ < 1. Besides 
the squeezing, the self-interaction 2kJ 2 also leads to the 
phase diffusion, which indicates a decay of the single- 
particle coherence. Such a kind of coherence is measured 
by off-diagonal elements of the single-particle density ma- 
trix — (a\aj)/N with i, j = 1, 2. Formally, one intro- 
duces the first-order temporal correlation function [l3T |: 



(AJ n 



1 2 

mm 
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(4) 



which is observable in experiments by extracting the vis- 
ibility of the Ramsey fringes @, ©, [M [10, El ■ One of the 
goals of this paper is thereby to present the relation be- 
tween the squeezing £ and the first-order coherence ■ 
Fet us first examine the exact numerical solutions of 
the time-dependent Schrodinger equation governed by 
the JFC Hamiltonian H2 . We consider that the spin sys- 
tem starts from the lowest eigenvector of J x , = 
e-™ J y/ 2 \j,-j), a particular CSS, Eq. ©, with 6 = tt/2 
and cf> = 7T. Such an experimentally realizable state can 
be prepared by applying a two-photon tt/2 pulse to the 
ground state with all the atoms in the internal 

state |1) [1, O EH. The spin state at arbitrary time t 
can be expanded as: |^) = J2 m Cm lii TO )i an d the am- 
plitudes c m obey ic m = e m c m + X_ m c m+1 + X m c m _i, 
where e m = 2nm 2 

X-4 



Cm(0) = (j,m\j 



. X m = + m)(j -m+ 1) with 

0. The amplitudes of the initial CSS are 

_ (_ 1) J+m t 2j ^1/2 



= ^j— (/ + >J ■ Obviously, 
C- m [U) = c m (0) for even N and c_ TO (0) = — c m (0) for 
odd TV, which gives the expectation value (J z (0)) = 
and the variance (J 2 (0)) = j/2. Note that some refer- 
ences adopt the Hamiltonian H3 = —VtJ x + 2kJ 2 to in- 
vestigate the BEC in a double- well potential [2(j, which 
corresponds to B.% for < case. In this work, we con- 
sider only positive O and k by assuming repulsive atom- 
atom interactions. If either f2 or k is negative, our results 
remain valid by using initial state \j,j) x [22| - 

Since c_ m (0) = ±c_ m (0) and X ±m = X Tm +i, we in- 
troduce linear combinations of the amplitudes Pm = 
c m ± c_ m . For even TV case, one can derive a closed 
m ^ ■ However, all Pm\o) = 
lead to p y m'(t) = 0, thus C- m (t) = c m {t). As a re- 
sult, dynamical evolution of the even N system is deter- 
mined solely by the equations of the amplitudes Pm with 
772 = 0,1, j. Similarly, for the odd N case, we obtain 



set of equations for p. 

„(-) 



III. QUANTUM DYNAMICS OF THE 
COHERENCE AND THE SQUEEZING 

The OAT model Hi (i.e., = 0) can be solved ex- 
actly in Heisenberg picture [l{, with its analytic results: 
A = 072)0' - 1/2) [l-cos 2 ^- 2 (4^)], B = -j(2j - 
1) sin(2«;t) cos 2j ~ 2 (2«;t), and C = j + A due to the time- 
independent variance (J 2 ) = j/2. The strongest (op- 
timal) squeezing £ = — (^S) 1 / 6 ^ -1 / 3 occurs 
at time t s ~ 6 1 / 6 A r_2 / 3 /2. The single-particle coher- 



ence #12 (i) = cos 
diffusion time Kid 



,N-1 



(2«t) 
(2/V)- 



-1/2 



-(*/*<*) with the phase- 
Obviously, ffi^td) = 



e _1 5j^(0) = 1/e. The coherence (t) has been mea- 
sured in experiment by extracting the visibility of the 
Ramsey fringe [12] . As shown in Fig. [lja), the opti- 
mal squeezing occurs within the coherence time due to 
t a < td- Moreover, the coherence decays to zero at 
to = 7t/(4k), and recover to unity at 2to [not shown in 
Fig.[T{a), see Refs.dd]. 

Except for N = 2, 3 cases, the JFC model (f2 7^ 0) can 
not be solved exactly [HI, [la, HI . Numerical simulations 
of the single-particle coherence g^ (t) and the squeezing 
^(t) are presented in Fig. [IJb)-(d) for N — 40 and vari- 
ous f2. Similar with the OAT case, the coherence g[^ (t) 
collapses to its local minimum at to then revives partially 
at about 2to, and the maximal squeezing occurs at t s for 
a small coupling £1 = 2k [Fig. [2(b)]. Two time scales to 
and t s tend to merge with the increase of Q. For f2 > Oo, 
both the coherence g^ and the squeezing £ reach local 
minima at the same time to [Fig.QJc) and (d)]. Here Flo is 
the optimal coupling to produce the strongest squeezing 
in the JFC model, such as fi = 4.2405k for N = 40. 

It should be mentioned that the spin state at to ex- 
hibits a very sharp probability distribution and a strong 
reduction of the number variance AJ Z [l(|. The loss of 
the coherence (or visibility) as an evidence of the number 
squeezing at time to, as shown in Fig.[TJb)-(d), has been 
observed by Orzel et al. @. Moreover, our results show 
that the phase diffusion is suppressed due to the appear- 



3 




between g^ and £ 



0.2 0.3 0.0 0.1 0.2 

t (k 1 ) t (k 1 ) 

FIG. 1: Time evolution of the squeezing parameter £ (solid 
thick lines) and the single-particle coherence g$ (t) (thin lines 
with open circles) for N = 40 and various Rabi frequencies: 
(a) fl = 0, (b) n = 2, (c) fi = n = 4.2405 (optimal coupling), 
(d) f2 = 8. Time i is in units of k -1 , and the units of Q is k. 



ance of number-squeezed state at £o ^ consistent with the 
experimental observations In experiments so far, 

however, the number squeezing were detected through 
the observation of increased phase fluctuation A<fi, or 
through an increased phase-diffusion time [24J. In fact, 
the number fluctuation has a nontrivial relation with the 
phase fluctuation, complex even for a single-mode light 
field [25l ]. Consequently, a more direct way to measure 
the number fluctuation A J z is necessary. 



IV. EXACT RELATION BETWEEN THE 
COHERENCE AND THE SQUEEZING AT THE 
MAXIMAL-SQUEEZING TIME 



e(t 



2<^(*o)) 



= 1 - 



J 



(6) 



which is valid for arbitrary f2. For instance = 0, 
g^2 = and £ = 1 at to = 7t/(4k); while for £1 > Nk, 

from Eq. ^ we get g^(to) — £(*o) — lj two trivial 
results due to weak squeezing. Hereafter, we focus on 
the coupling around its optimal value Oq, with which 
the strongest squeezing £ = £(i ,f2o) can be obtained 
at time to- According to Eq. ([5]), one can measure £o by 
readout the coherence g[^ in atomic interference experi- 
ments BSIMClIClI. 




FIG. 2: (color online) Husimi Q function in the phase space 
(4>,p^,) for various times: (a) t = 0, (b) t = 0.045, (c) 
t = to = 0.09064 (the optimal squeezing), (d) t — 0.135, 
(e) t = 0.183127, (f) t = 0.225, (g) t = 0.26988, and (h) 
t — 0.315, corresponding to the times indicated by the ar- 
rows in Fig.fTJc), respectively. Other parameters are N — 40, 
Q, = f2o = 4.2405k, and time t is in units of kT 1 . 



There exists an exact relation between the coher- 
ence and the squeezing £ at time t . To see 
this, let us examine the Heisenberg equations: J x — 
—2K(J z Jy + J y J z ) and J z = flJ y . The first equation 

gives the relation between the coherence g^ 2 and the 
squeezing angle 9 min 

j t g[l } = 2«Aj- 1 tan(2# min ). (5) 

As g-^ reaches its local minimum at to, (dg^ /dt) to = 0, 
then the squeezing angle m ; n = (or n) provided A^=Q. 
Combining the two Heisenberg equations, we obtain fur- 
ther dJl/dt = Sl(J z J y + J y J z ) = -{n/2n)dJ x /dt, which 
yields (J z ) = A — Q(J X ) /(2k) with the integral constant 
A. For the initial CSS \j, —j) x , we have A = — Q/k)/2. 
At time to, B — and A > 0, and the minimal variance 

( AJ n)Ln = <£ ~ A )/ 2 = ( J z) i Le -> the squeezing along 
the z axis [17| . As a result, we obtain a simple relation 



The time scale to can be obtained based upon the 
phase model [l6|. By replacing J z — > p^ = —id^ and 
Jx —> j cos <j), the JLC Hamiltonian H2 can be rewritten 
as = 2kp^ + flj cos (f> [IH , where <j) is relative phase 
between two bosonic modes. The phase model allows us 
to treat the spin system as a simple classical pendulum 
oscillating around the minimum of the Mathieu potential 
cos <f>, i.e., cf> — tt. In the Josephson regime 1 < Q/k -C N, 
the pendulum rotates in the phase space (4>,p$) with the 
effective frequency u) e g — \/2kQN [H, H3] ■ To illustrate 
this motion, we calculate the Husimi Q function 

Q(0,$ = KM*(*)>| 2 (7) 

in the phase space (4>,p$), where p^ = — jcos((9) de- 
scribes the population imbalance of the two modes, and 
the polar angle <j> represents the relative phase [28f . The 
CSS \9, 0) is given in Eq. ©. As shown in Fig. Ufa), the 
Q function is a circle for the initial CSS, which represents 
Poisson distribution of the number variance (AJ Z ) 2 and 
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the phase uncertainty (A0) 2 [12j| , As time increases, it 
becomes an elliptic shape [Fig. [DJb)], rotating clockwise 
in the phase space. After a duration to, the ellipse elon- 
gates horizontally corresponding to the optimal squeez- 
ing of (AJ Z ) 2 [Fig.[2fc)]. It seems reasonable to suppose 
that the motion of the ellipse is consistent with that of 
the pendulum. In fact, the trajectory of the pendulum is 
just passing through the horizontal axis (p^, = 0) at time 
T/4, where T = 27r/w c ff is the period of the pendulum. 
As a result, we get the maximal-squeezing time [l6[ 

K * ^ I = _y_ (8 ) 

which is valid for large N (> 10 3 ). As shown in Fig.^e)- 
(h), for t > 2to (~ T/2) the Q functions almost recover 
to original shapes, due to partial revival of the squeezing 
£ and the coherence [see Fig. QJc)]. 




N 

FIG. 3: (color online) Numerical simulation of the optimal 
coupling Qo as a function of N (solid circles) in the normal 
scale. The solid line is the fit to power law f2o = aN b with 
a = 1.07827 and b = 0.32655. The inset: comparison linear fit 
(blue line) of numerical results (open circles) with Qo = N 1 ^ 3 
(red line) in the log-log scale, fio is in units of k. 



V. OPTIMAL COUPLING AND THE 
STRONGEST SQUEEZING 

To create the strongest reduction of (AJ Z ) 2 as shown 
in Fig. dfc), we need to determine the optimal coupling 
Slo as a function of particle number N. Note that the 
optimal squeezing occurs at t s for the OAT and to for 
the JLC, respectively. For large N, the latter time scale 
should be comparable with the former one as £1 = fio- 
Such a non-rigorous comparison enables us to suppose a 
power law as Qo/n ~ A 1 / 3 . Numerical solution of ilo is 
presented in Fig.[3]for N up to 2 x 10 5 . We fit the data as 
S1 /k = aN b and find the power-exponent b = 0.32655, 
very close to the expected value 1/3. From the inset of 
Fig. we also find that the larger number N is adopted, 
the better fit is obtained. 

In Fig. HI we investigate the optimal squeezing £o as a 
function of N. The fitting result is £ - 0.85787V- 1 / 3 , 
which is slightly smaller than the OAT result £o ~ 
(4/3) 1 / 6 A~ 1 / 3 ~ 1.0491 A" 1 / 3 . Small difference of £ be- 
tween the OAT and the JLC does not deteriorate the ad- 
vantages of the latter scheme. In fact, there is no number 
squeezing in the OAT model due to [Jf , Hi] = 0. In our 
case, one can realize (J^(t )) = £q(^z(0)) wrtn £o < 1, 
indicating the appearance of the number-squeezed state 
Such a kind of squeezed state has been observed 
in optical lattices @, optical trap and atom chip 
However, the observed squeezing £o ~ 0.1 for 
N = 4x 10 5 weaker than our result £ - 1.467 x 10" 2 
for f2 = 58.05k and N = 2 x 10 5 . Finally, within inset 
of Fig. [4} we show the time scale to as a function of N. 
Inserting VL /k = A 1 / 3 into Eq. ©, we find that analytic 
expression of to gives good agreement with the exact nu- 
merical simulations. 
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FIG. 4: (color online) The optimal squeezing £o for the OAT 
(open circles) and the JLC (solid circles) as a function of N. 
The blue line is given by £ ^ (4/3) 1/6 A^- 1/3 ~ 1.04917V~ 1/3 ; 
the red line is a fitting curve as £o = 0.8578iV" 1/3 . The inset: 
numerical simulations of to (open circles), and the analytic 
result of Eq. ([8| with Qo/k = N 1 ^ 3 (solid line) as a function 
of N for optimal squeezing. The time to is in units of k," 1 . 



VI. CONCLUSION 

In summary, we have investigated optimal spin squeez- 
ing in a two-component Bose-Einstein condensate with a 
Josephson coupling. We show that: (i) the squeezing £ at 
time to aligns along the z axis, which is equivalent with 
the number squeezing d, [^, [lfl, HH, fl7|, an d is desirable 
for high-precision atom interferometry [28|, |2!| ; (ii) there 
exists a simple relation between the squeezing £ and the 
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single-particle coherence g[^ at to, Eq.© and Eq.©, 
from which it is possible to measure the number vari- 
ance A J z by readout the coherence g[^ in the interfer- 
ence experiments; (iii) the strongest squeezing with the 
power law £o — 0.8578A^ 1 / 3 is achievable by applying 
the optimal coupling VIq/k, ~ N 1 / 3 . We also discuss the 
maximal-squeezing time to via the phase model and the 
Husimi Q function, and find that analytic result, Eq. l[5]). 
agrees with its numerical simulations. 
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